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O ' I The problem of color superconductivity in dense QCD is reconsidered in the improved rainbow 

approximation to the Schwinger-Dyson equation. The effect of the unscreened magnetic modes of 
gluons on the value of the color condensate is studied. In particular, it is shown that, at sufficiently 
large values of the chemical potential, these modes lead to the enhancement of the superconducting 

^ ■ order parameter. The interplay between the instanton induced interaction and the one-gluon induced 

O ' one in color superconductivity is discussed. 
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I. INTRODUCTION 



> 

00 

^T I Quantum chroniodynamics (QCD) has been remarkably successful in describing the interactions of quarks and 

gluons at short distance scales. In this distance region, QCD is a perturbative theory due to asymptotic freedom [nl. 
At large distances, on the other hand, the effective running coupling of QCD becomes strong and non-perturbative 
approaches, like sum rule methods H, lattice computer simulations 0], instanton computations B, and Schwinger- 
Dyson equations for Green's functions ||^ should be utilized to analyze the theory. 

There exist, however, some extreme conditions where QCD is yet to be tested. An example of such an extreme 

condition is cold quark matter at high density. There have been attempts to investigate dense QCD using perturbation 

O i' theory [gj, instanton calculations 0-^ and lattice simulations |1^. Such dense matter may exist in the interior of 

1^ I neutron or the so-called strange stars UM, with baryon number densities exceeding a few times the normal nuclear 

'^ ■ density uq ~ 0.17 fm~ . Besides that, the dense enough quark matter could be created in accelerators by heavy ion 

^ , collisions. Therefore, the study of such a system is not of pure academic interest. 

At high density, it is believed that the quarks form a Fermi surface in a very similar way as the electrons do it in 
^ . , metals. Following this similarity further, it is natural to ask whether there is an analogue of the superconductivity 
j^ ' in the cold quark-gluon plasma. As in nretals, the presence of the Fermi surface should considerably increase the 
density of states of low-energy quasiparticles. As a result, an arbitrarily small attractive interaction in the diquark 
channel would be sufficient to create the Cooper pairs and, thus, to produce the (color) superconductivity. In Ref. [||, 
it was shown that an attractive interaction in the color-antitriplet diquark channel indeed appears. As a result, a 
non-trivial superconducting order parameter develops and the color gauge symmetry SU(3)c gets broken (by the 
Higgs mechanism) down to its SU{2) subgroup (in the case of two light quark flavors). 

Recently, interest in the study of the color superconducting phase considerably increased ]7|-P,p^-p2[ . The renewed 
interest was triggered by the observation that the necessary attractive interaction could be exclusively due to the 
instanton effects iR-Pl- In addition, some other fascinating features of the broken phase, such as color-flavor locking and 
a new type of chiral symmetry breaking p2| , were revealed. From a technical point of view, it was also demonstrated 
that the renormalization group method is an alternative and a very efficient tool in studying the general properties 
of dense QCD |l|-|l|. 



* Permanent address. 
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While it was believed that the color superconductivity was quite well understood at least qualitatively, the recent 
results of Refs. |1^,Q seem to indicate that there is an inconsistency in all the previous approaches, based on 
the straightforward use of the simplest Nambu-Jona-Lasinio model and the Bardeen-Cooper-Schrieffer (BCS) type 
of analysis. The key observation in Refs. |ljj2y] is that the long range interaction, mediated by the unscreened 
gluon modes of the magnetic type, may considerably enhance the value of the superconducting order parameter. 
Such an enhancement would be of great importance in studies of neutron star properties. Note that a somewhat 
similar conclusion about enhancement was also reached in Ref. pl[ | where the effective low-energy action of the quark 
quasiparticles around the Fermi surface was derived. 

Being motivated by that observation, in this paper, we reanalyze the problem by using the conventional Schwinger- 
Dyson (SD) approach. More precisely, we study the SD equations in dense QCD in the so-called improved rainbow 
approximation in which the full vertices coincide with the bare ones. Such an approximation includes the one-loop 
polarization effects in the gluon propagator and, thus, is the simplest non-trivial approximation that takes into account 
the screening effects. Our expression for the fermion gap derived in the SD approach turns out to be essentially the 



same as in Ref. 1 20 where the renormalization group method was mainly used. 



This paper is organized as follows. In Sec. ||, we present the general structure of the quark propagator that 
accommodates the color superconductivity. In Sec. |III| , we discuss the structure of the gluon propagator with the 
one-loop polarization effects taken into account. The role of the Meissner effect in the color superconducting phase of 
QCD is discussed in Sec. IV. Then, in Sec. M, we derive the SD equation (in the improved rainbow approximation) 



for the quark propagator and present an analytical estimate for the solution. The discussion of the results is given in 
the concluding Section VI. In Appendix H, some formulae used in the derivation of the gap equation are considered. 



And, in Appendix M, we present an approximate analytical solution to the SD equation. 

II. THE MODEL AND NOTATION 

Here we consider QCD with two light quark flavors {u and d) in the fundamental representation of the SU{3)c color 
gauge group. At sufficiently large values of the chemical potential, the current masses of the quarks can be neglectedfj. 
As a consequence, the model is invariant under the (global) chiral SU{2)l x SU{2)r transformations. 

Following the findings of Ref. |^, we are interested in studying the diquark condensates e"^^ £ab-i{{i^a)'^ Cij^l) ^i^d 
e'^^ eabsiii^l)'^ 0^51^1) , responsible for the color superconductivity (here a, b and i, j are the color and flavor indices, 
respectively). We note, that in perturbation theory the two condensates are related by Ua{^) symmetry. Nonpertur- 
batively, however, the instanton effects are expected to break the latter symmetry [RH]. 

In the presence of a large chemical potential, it is commonly assumed that the quarks near the Fermi surface (or 
rather low-energy quasiparticles) are weakly coupled due to asymptotic freedom [g^. Then, it is natural to analyze 
the dynamics of diquark pairing by using the SD equations with a perturbative kerneln pl| . At this point, it is 
appropriate to mention that the realistic value of the chemical potential for the quark matter in the interior of a 
neutron star is expected to be somewhere in between 200 and 700 MeV (so that the baryon number density is between 
about no and lOno). Then, the corresponding value of the coupling constant as should be between about 0.4 and 0.8. 
At these values of a^, our calculations, based on the perturbative kernel, may just start to fail. At the same time, 
since the values of the coupling constant are not too large, it is plausible that the orders of the magnitudes of the 



main results could still be trusted. We will return to the discussion of this point in more detail in Sec. VI. 

Here we should mention that the SD equation approach to the study of color superconductivity in dense QCD 
was also used in Refs. pjl^,^,p3. In this paper, however, we pay special attention to the effect of the unscreened 
magnetic modes of gluons. 

Instead of working with the standard four component Dirac spinors, in our analysis below, it is convenient to 
introduce the following eight component Majorana spinors: 



^/2 V^' 



^As is argued in Ref. |l2], the strange quark is also sufficiently light to be included in the model. However, we are not interested 
here in the specific effects of the third fiavor. Our prime goal here is to clarify the role of the long range gluon interaction and, 
thus, it is sufficient for our purposes to deal with the simplest model of two light flavors. 

■^Note that according to Ref. flSl, the effect of instantons could be negligible compared to the one-gluon exchange at sufficiently 
large chemical potential. 



where C is a charge conjugation matrix, defined by C~^7^C = — 7J and C — — C^. In the new notation, the inverse 
fermion propagator (defined over the true vacuum) that accommodates the possibility of the diquark condensation 
reads 



(2) 



^^^'' \ A A{p)[{po-fih' + B{p)flJ' 

where A — •^'^A^j'-' and A^-'j, = e''^eab3^3 (note that A3 is a matrix in the Dirac space). By definition, f= —p- 7. 
Now, after inverting the expression in Eq. (0), we arrive at the following propagator: 

G(r,)=^( ^iW^' -[{po + f^h" + B{p)f]-'A{p)-^AR2{p)-'\ 

\-[ipo-^lh" + B{p)^]-'Aip)-^AR^ipr^ R^ip)-^ J' 

where 

i?i(p) = A{p) [(po + A^)7° + Sb¥] - A [{po - /i)7" + B{p)fi] "' A{p)-'A, (4a) 

R^ip) = A{p) [ipo - /i)7° + Bip)f] - A [{po + m)7° + B{p)f] ~' A{p)-^A. (4b) 

In Sec. |y|, we use the expression for the propagator in Eq. (y|) as an anzatz for the solution to the SD equation. 
Notice that, with the given choice of the condensate, the most general structure of the wave function rcnormalizations 
A{p) and B{p) (in flavor and color spaces) has to be as follows: ^;lfy{p) — a[p){5ab ~ SazS^^jS^^ + a{p)5a3,5b3,5'^^ and 
^abip) ~ b{p){Sab ~ Sa35b3)S^^ + b{p)Sa3Sb3S'^^ , whcrc a(ji) , a{p) and h{p), b{p) are some scalar functions. 

III. VACUUM POLARIZATION IN THE GLUON PROPAGATOR 

In this section we discuss the vacuum polarization effects in the gluon propagator. 

The calculation of the polarization tensor was performed by others ||25[-p7t and we are not going to repeat it here. 
Instead we present the final result and comment on its essential features. 

Let us start from the expression for the inverse gluon propagator in a covariant gauge with polarization effects 
taken into account. Its explicit form is given by 

[V^^{k,, k)) "' = ^5^''k^P^, + ^5^''^Pl + ^<5^^^,.(fco, k). (5) 

where the color indices A,B = 1, 2, . . . , 8, and the standard projection operators P^ and PJliy, 

P^.^=9^.u~-^, and Pjl^^-^, (6) 

were introduced. The gauge fixing parameter d in Eq. (^ is arbitrary. 

At a finite density of quarks, the Lorentz symmetry is explicitly broken and, as a result, the polarization tensor 
H^u could also contain the third tensor structure [ p8| : 



fc2 {u ■ k) {u- ky 






where m^ = (1,0,0,0). This last tensor, similarly to P^, is transverse, k^PJ^^^ = 0. The set of all the three tensors, 
introduced so far, satisfy the following set of multiplication rules: 

p±p±^p±^ Pllpll=pll, p»p» = ( ^ 1 )P", (8a) 

\{u-kY J 

p±p\\ ^ plIp-L ^ Q^ pWpu ^ pupW ^ 0, P^'P^ = P^P"" = P", (8b) 

where we assume that the operator products are defined by the appropriate contractions involving the Minkowski 
metric. While working with the gluon propagator, it is much more convenient to use the following set of mutually 
orthogonal projection operators: 



0(1) = pJ 



(u- ky 



o(i)o(i) = o(i) 



P", 0^2) = - 



(u-fc)2 



-p«, (9(3) ^ pll 



0(1) + (9(2) + (9(3) ^ 



fiiy 



' ^.y 



fiiy 



9^^i 



0(2)0(2)^0(2), 0(3)0(3)=0(3), 

O(0o(j)=0, for i^j. 



(9a) 

(9b) 
(9c) 



At finite chemical potential, in the one-loop approximation, the calculation of the polarization tensor reduces to 
evaluating the contribution of three diagrams: with gluons, quarks and ghosts running in the loop pq |, respectively. 
In the hard dense loop approximation, the result for the polarization tensor is given by |25|-|27J| , 



n"0(fco,fc) = ni(fco,fc), 
n"*(fco,fc)-fco^ni(fco,fc), 

|/c|2 

ff^(fco,fc) ^ U^^- - ^ ) n,(fco,fc) + ^^Mko,k), 



\k\' 



where 



ni{ko,k) = 2M^{ -^In 
nt{ko,k)=M^ 



kg 

m 



\k\ 



m- 



kg-\k\ 
Mko,k). 



\k? W 



fco 



- 1 - i7r^^6i(-fc2) , 
2|fc| V 



(10a) 
(10b) 

(10c) 



(11) 
(12) 



Here we use the notation NP = {gsjJLyjNf /2n)'^ . What is remarkable about this result is that it coincides with the 
polarization tensor derived in the framework of classical transport theory of dense Yang-Mills plasma pQjp^ • 

Notice that the polarization tensor in Eq. ( |l0| ) has a nonzero imaginary part for space-like gluon momenta. This 
imaginary part is responsible for the so-called Landau dampingof the gluon field with space-like momenta. Also, it 
is responsible for the quark damping around the Fermi surface [E5 29 . In our analysis below, however, we neglect the 
effects of quark damping. Such an approximation is partly justified by the fact that the damping rate goes to zero 
linearly as one approaches the Fermi surface |29|. 

The gluon polarization tensor in Eq. (|lO|) is transverse 



fc^n^,(fco,fc) = o. 



(13) 



Because of this last property, it is natural to try to express the real part of the polarization tensor in terms of the two 
transverse projection operators, O^^^ and O'^', introduced earlier. This turns out to be an easy task and the result 
takes the following nice form: 



Hp. = -o^l^n, + 2o(,i) (n* - M^) . 

With this expression at hand, we rewrite the inverse gluon propagator in Eq. (0) as 
{V^^^ik,, k)) ' = ^5^^ {e - no OW + ^5^^ {^ + 2n, - 2M^) O^^^ - 



-^^^^o(i). 



(14) 



(15) 



Then, by making use of the properties of the projection operators, we could easily invert this expression and arrive 
at the final form of the gluon propagator that includes the (one-loop) screening effects 



2^;:.(fco,fc) = 



xAB 



1 



fc2 



_o(i) 

Tit '"' 



iS^"" 



1 



fc2 + 2Ilt 



0(^)-z(5^^ — 0(3) 



(16) 



This representation of the gluon propagator, in addition to its convenience for calculations, also has another advantage. 
It allows to separate different gluon modes in a very simple way. Indeed, to make such a separation, we just need 
to project the general gluon state by using one of the three operators, 0^^\ O(^) or 0^^\ respectively. As a result 



S^) 



O^^^A" (magnetic), a^'^ 



we come to the following types of modes: a/i 

(longitudinal). The name "magnetic" here reflects the point that the corresponding projection operator does not have 



O^ii^A" (electric) and a^}^ = O^^} A" 



the electric components OIq — OqJ = 0. Apparently, the rest of the transverse modes are the "electric" modes. As 
for the unphysical longitudinal mode, its definition is standard. 

In the rest of this Section, we will consider the part of the gluon propagator connected with the unbroken, SU{2), 
subgroup of the SU{3)c, i.e., A,B = 1, 2, 3. The part of the gluon propagator connected with the broken generators, 
A, B = 4, ... ,8, will be considered in the next Section. 

Because of a specific dynamics of diquark pairing in the vicinity of the Fermi surface, it is believed that the most 
relevant gluons mediating the interaction are those with the space like momenta. Indeed, when the quarks around 
the Fermi surface scatter, their energy does not change much, while their momentum could change for as much as 
2pp. Thus, we could mimic the screening effects by the following asymptotic form of lit in the region |A;| ^ fcp: 

nt(fco, k) c 2M'S- (l-lS + ..]+ znM'X (l - S-] ei\k\' - kl). (17) 

|fc|2 \^ 3|A:|2 J 2\k\ \ \k\y 

In fact, this is a very good approximation even for the values of the ratio fco/|A:| as large as 0.8 at which the deviation 
of the asymptote from the exact expression is of order of 5%. If we keep only the leading term in Eq. (|l7|), the same 
magnitude of deviation is reached at about fco/|fc| ~ 0.4. 

By substituting this asymptote into the gluon propagator, we see that the magnetic mode aj/" produces the 
long-range interaction, while the electric mode a)^ gets screened out. In what follows, it is convenient to use the 
representation of the gluon propagator in terms of spectral densities: 

1 f°° dz 
m — / 

-0 2tt J_^ ko- z + ie 



1 /"CXD ^ r 

iV^^^iko, k) 0, S^^ hm - / '—— p^{z, k)0^^) + p,(z, k)0^;^) + pii (z, k)0^^^ , (18) 



where 



pAz. k) ^ ^^^l!flLe(|fcp _ ,^)^ (19) 

Pe{z, k) ~ 27r sgn(z),5(z2 - jfcp _ 21/^), (20) 

p\\{z, k) = 27rd sgn{z)5{z'^ ~ \k\^). (21) 



Therefore, by making use of the representation in Eq. (Ilq), we arrive at the following expression for the gluon 
propagator in Euclidean space (fco — ik/C): 

""^^ ' \k\^+7rA'P\ki\/2 """ kl + \k\^ + 2M^ ^" fc| + |fc|2 ''''' ^ ' 

where A^B = 1,2,3. This propagator correctly describes the gluons in the one-loop approximation in the soft 
momentum region. In Sec. M, we shall use it in the SD equation for the quark propagator. As we shall see there, 
because of the long range interaction mediated by the first term in the propagator in Eq. ( p2| ) , the result for the order 
parameter is going to be quite different from that obtained in the theories with local interactions [ [l2||l8| . Notice, that 
according to the arguments of Refs. |19,M, the magnetic modes of the unbroken SU{2) subgroup of SU{3)c should 
develop no screening even after taking into account nonperturbative effects. 

IV. THE MEISSNER EFFECT 

In this section, we study the role of the Meissner effect in color superconducting phase of dense QCD. Based on 
pure symmetry arguments, it is clear that when the color S'f/(3) symmetry spontaneously breaks down to SU{2), five 
out of eight gluons should get masses by the Higgs mechanism. 

To estimate the value of the Higgs-like gluon mass, we need to calculate the appropriate contribution to the vacuum 
polarization tensor that results from the non-diagonal term in the fermion propagator in Eq. (pf). To get a rough 
estimate of the mass, we could completely neglect the effects of the wave function renormalizations (in the next 
Section, we will give a justification for that). Then, we get the following order of magnitude expression: 

V(0) ^M^o,^ f d^Qdqifi^\A^-^\^ ^ a, , 



where the order parameter A'^~) is defined in Sec. ^ [see Eq. (p2[)]. This seems to suggest that the five gluons obtain 
relatively large masses (of order asfi^) due to the Meissner effect. As we show in a moment, however, this conclusion 
is not quite right. 

Before going into more details, it is instructive to remind a few facts about the Meissner effect of the low-temperature 
superconductivity in ordinary (non-relativistic) metals |3C| ]. As is well known, there exist two characteristic scales in 
the theory: the coherence length, ^, and the London magnetic penetration depth, A^. For our purposes here, it is 
sufficient to recall that while the coherence length, as T ^ 0, is directly related to the value of the superconducting 
order parameter (^ ~ 1/|A|), the London penetration depth is independent of A and is given in terms of the mass 
and the density of the electrons (A^ ~ yjrnjn^i). 

The London penetration depth Az, coincides with the actual magnetic penetration depth A only in type II su- 
perconductors with Ai, ^ ^: I/Al coincides with the value of the running Higgs-likc mass of a plasmon at zero 
momentum and the region of small momenta yields the dominant contribution to the actual penetration depth in 
type II superconductors. On the other hand, in type I superconductors, it is the running Higgs-like mass at momenta 
A <C fc <C /^ that yields the dominant contribution to A. In this case A = \p where Ap is the Pippard penetration 
depth0: Ap ^ {Xlifl^ @,|I|. 

Regarding the cold dense quark matter, it can hardly be a type II superconductor because, as we shall see from 
our final estimate for the gap in Eq. (p3), |A(~'| is |A*^^^| ^ Mq, implying that Al ^ ^. Therefore, most likely, it 
is a type I superconductor. This means that the actual penetration depth of a magnetic field is given not by the 
London expression but by the Pippard one. Technically the difference comes from the fact that the zero-momentum 
expression in Eq. (E3h , related to the definition of the London penetration depth, is valid only in a very small region 
of momenta, |fc| <C lA^^^^I. In the most important intermediate region, jA*^^^! ^ |fc| <C /i, however, the expression has 
the following asymptotic behavior |pl|: 

V(k)^Ml^^. (24) 

|fc| 

Since the relevant momenta for the Meissner effect are of order A~^ 3> |A^~)|, from Eq. (p4), we indeed see that the 
penetration depth is equal to Ap ~ (Mg |A(~)|)~^/'^ rather than A^ ^ l/Mg. 

Now, what the role does the Meissner effect play in the perturbative kernel of the SD equation? To answer this 
question, let us consider how the corresponding propagators of magnetic gluons are going to be modified by the 
Meissner effect. Apparently, we could still use the representation for the propagators as in Eq. (ttq). The spectral 
density, though, should be modified appropriately: 

pt'\^S)^^ ^^^^^^-^ m?-z\ (25) 

(|fc|3 + Af2|A(-)|)2 + (^/2)2M4z2 

where the superscript M stays here for the "Meissner effect". After simple calculations, we arrive at the following 
expression for the corresponding propagator in Euclidean space (fco — ik/^)'- 

^Vlf^^^i^h,k)^-6^^ '^1 ^^0W + ..., for A, i? = 4, . . . , 8, (26) 



where the ellipsis denote the same electric and longitudinal contributions as in Eq. (E2[) . By comparing the propagators 
of magnetic modes in Eqs. (g^ and (|26|), we see that the Meissner effect could be accounted by the formal replacement 
^4 -^ fc4-f c|A^~)| with c = 0(1) in the magnetic term of the five propagators in Eq. (Eq). Therefore, it is obvious that 
the difference between the propagators in Eqs. (E2h and (Eq) could become important only in the range of momenta 
|fc4| ;$ lA^^^^I and \k\ ^ (Mq |A'^~^|)^/'^. In the next section, we shall see that the mentioned region of gluon momenta 
is not large enough to modify the leading asymptote of the solution to the gap equation. 



^We would like to thank S. Hsu for pointing our attention to the difference between the definitions of the London and the 
Pippard penetration depths. 



V. SCHWINGER-DYSON EQUATION IN DENSE QCD 



In this section we derive the SD equation in dense QCD, ignoring the masses of quarks. AU the necessary constituents 
of such an equation were given in the previous three sections. In the improved rainbow approximation, the equation 
reads 



(Gabip))-' =(G^at\p) 



+Anas 



Aira^ 



dS fj 



(27r)4 l^ -7^ 



/ , / , T^'aGa'b'{q)Tf,>b 



a'b' A=l 



r 







-7^^ 



V 'Df.uiq-p) 



(27r)4 V -7^ 



A^ )EE^a^^G,,,,(9)r,-f 



a'b'A=4 



Y 



\[^ _^.W^\<i~pi 



(27) 



where V^^^k) is the propagator of gluons that correspond to the unbroken SU(2) subgroup [see Eq. (|22|)], while 

'Dliv {k) is the propagator of those five gluons whose magnetic modes are modified by the Meissner effect [see Eq. (p6|)] . 
Notice, that the overall factor 6^^ is omitted in the definition of the propagators here. Regarding the rest of notation, 
Gab{p) is the full fermion propagator, and G)^^{p) is the perturbative one. With our choice of the order parameter 
orientation in the color space, A^b ~ Cabs , the explicit form of the generators of the unbroken color subgroup reads 



T' 





I , where A= 1,2,3, 



(28) 







and a are the Pauli matrices. Then, by making use of the well known identity for the Pauli matrices, we arrive at 
the following result for the summation over A: 



3 1 1 

E "^^oT^b = T: i^a'b - Sa'sSbs) {Sab' - SaS^b's) ~ T ^^'^'^' ^ '^aS^a's) {hb' - ^bS^b's) 



A=l 

By making use of the identity for all eight generators of SU{3): 

8 



Err^A rp 



A 
b'b 



1 



:Sa'bSab' 



:Saa'Sbb', 



(29) 



(30) 



A=l 



we also obtain the following one: 



* 1 1 1 1 

22 "^^a^Wb = -r::Saa'5bb' + 2 [Sa'bSazh'Ti + 5 ab- 5 a- Z^bS.) - T {5 aa' 5b?,5b' Z + hb'Sa^Sa'z) - T'^aS'^M'Ja'S'Jb'S, (31) 



A=4 



Now, upon inserting the identities ( |29| ) and (31) along with the expression for the fermion propagator (^) into 
Eq. (p7|), we arrive at the set of integral equations: 

ab)A(-n.)^^/^^(9P..(.~p)-eH.-.); 



(2^)4 a{q) 
A'^-\q) Tr (7^A^+^7'^A^^^ 



A(+)(g) Tr (7^A^~^7''A^^^ 



g2-(%1-/i)2-|A(-)|2 qi-{b\q]+^ir-\Ai+)\^ 

(a-l)ipo + f^)Tiab-l)\p\^~^J^^^[3V,4q-p)+Vl'JHq-p) 



(32) 



(90 + b\q] - fi) Tr (7"7^a[+^07-a(±)) (go - b\q] - ^I) Tr (707^A(-)7%''A(^^) 



{a - l){po + fi) T {ab - l)\p\ 



g2-(fo|9l-/i)2-|A(-)|2 
27ra., f d'^q 1 



gg-(6|gH-/i)2-|A(+)|2 



(33) 



3 J (27r)4a(g) 
Tr (707^A^+)707'^A^^^ ) Tr 



-D^^JHq-p) 



qo -b\q\+ fi 



qa + b\q\+ ^j. 



(34) 



Here, Tr denotes the trace in Dirac indices, and A^^^ are complex functions of momentum, defined by either^ A^lp) - 
a{p) A(+)(p)A^+' + A(-)(p)A^"^ (parity odd order parameter) or A^ip) = 0(^)75 A(+)(p)A^+^ + A(-)(p)A^"^ 
(parity even order parameter), while 

are the free quark (antiquark) on-shell projectors ||l9| , |3^ , |33| p|. 

Let us start our analysis by considering the equations for the wave function renormalizations ( p3| ) and (p3) first. In 
order to get a rough estimate for the wave function renormalizations, it is sufficient to substitute a = b — a — b— 1 



along with A*^^^ = into the right hand side of Eqs. ( p3| ) and (34). Then, it is straightforward to show that (i) the 
wave function renormalization of the temporal and the spatial parts of the fermion kinetic terms are not equal, and 
(ii) the one-loop corrections to a{p) and b{p) develop logarithmic divergences at p4 = and \p\ — p,. As is clear, 
these logarithmic divergences at the Fermi surface are going to be removed when a non-zero order parameter A^"' is 
reintroduced. Notice that, because of the Meissner effect, no infrared divergences develop in a{p) and b{p) either. At 
the end, we obtain the following estimates: 

a,b,a,b '^ 1 + Const • ag (/i) In ,_, . (36) 

By taking into account the expected value of the order parameter [see Eq. (E2h below] , we actually see that all wave 
function renormalizations are close to 1 if the coupling constant a.s{pi) is weak. Therefore, we conclude that, in the 
leading order approximation, it is justified to neglect the wave function renormalization effects. 

The fact that a«6~a~6«l considerably simplifies the study of the gap equation (£^). However, one still 
has a rather complicated set of two coupled integral equations for A^") and A(+', respectively. At this point, it is 
important to notice that only A'^^ defines the gap in the quasiparticle spectrum around the Fermi surface. Also, in 
the leading order approximation, it is sufficient to keep only those terms on the right hand side of Eq. (p3) which 
become singular at the Fermi surface as A^^' -^ 0. Then, we see that the equation for A^~\p) decouples: 

irn f d^a / , A(-)((7)Trf7^A^+^7'^A^"') 

We further simplify the gap equation by assuming that A^'{p) = A^ ' {pi^pp) with pp — (0,0,//). This approx- 
imation with the order parameter being a function of only p4 is partly justified by the structure of the magnetic 
(dominant) part of the perturbative kernel in the SD equation. Indeed, in the vicinity of the Fermi surface, the 
dependence of the kernel is more sensitive to changes in p^ than to changes in \p\ . This is due to the fact that while 
the dependence on p4 comes from the linear term in the denominator of the propagator of the magnetic modes in 
Eqs. ( p^ ) and ([2q), the dependence on \p\ comes from the cubic term. 

In order to linearize Eq. (|37|), we substitute A'^^\p) — > A'~) = A'^~\p4)\p^^Q in the denominator. Then, after 
taking the trace over the Dirac indices and switching to the Euclidean space (go — *'?4)i the integral over the angular 
coordinates in Eq. ( p7[ ) can be done exactly (see Appendix |AJ for details). After performing the subsequent approximate 
integration over the absolute value of the spatial momentum, we finally arrive at 

A("-'(n4 =^ / dq^ ^^^^ Qln -L±2 ^ _ In ^, ,, ,, ,, -|- 12 n ^-^ -f 12d 

(38) 

As is easy to trace back, the first two terms in this expression come from the interaction mediated by the magnetic 
modes: while the first one is connected with the three gluons of the unbroken SU{2) subgroup, the second term 
comes from the interaction mediated by the five gluons subject to the Meissner effect. The third term in Eq. ( |3q ) 
is connected with the electric modes, and the last term comes from the longitudinal modes. By noticing that the 



^Notice that we factored out the renormalization of the wave function a{p) in the definition of A3(p). 
^In our notations, A'~' is the same as (Ai)* of Ref. |33]. 



main contribution to the right hand side comes from the region of momenta jA^ )| <C Iqa] <IC jU, it is obvious that the 
Meissner effect is of no importance in this gap equation. Therefore, by putting Mq = 0, we arrive at 

where A = e^'^/^2^f/{V2T:M'^) = e^''/^AT:)^/^fi/a^/^. 

Notice that, because of the absence of the logarithmic factor in front of the gauge parameter in Eq.(B8|), the 
longitudinal giuon modes become relevant only in the next-to-leading order. Here, however, we keep their contribution 
in order to estimate the size of possible corrections to the magnitude of the order parameter due to subleading effects 
[see the discussion following Eq. (|4^)]. The contribution of the electric modes is also subleading. However, because of 
the large logarithm ln(/z^/M^) ~ ln(l/Q;s), their contribution to the gap is important at asymptotically high density 
of quark matter when the couphng as(^) is weak [see the discussion after Eq. ([i^)]. 

To get a rough estimate of the order of the gap, we could just put a sharp ultraviolet cut-off at (74 = /it and then, 
assuming that A^^'{q4) sa A^^' — Const in this region of energies, we come to the algebraic equation: 

ic^^r ^In^.^ln^^, (40) 

leading to the following estimate for the gap: 

|A(-)|^(l^e-/.expr-3,/^V (41) 



a 



5/2 



A more rigorous solution of the integral equation ( p9| ) is presented in Appendix H. The latter gives a slightly smaller 
value of the gap: 

where e = 2.718 . . .. Notice that the exponent in this expression is completely determined by the magnetic modes. 
The main role of the electric modes is in changing the power of Qs in the prefactor: from gj'^ to gj^- At last, the 
contribution of the longitudinal modes in the prefactor is formally of order one and subleading, though it is essentially 
different from 1 for large values of the gauge parameter jdj. Obviously, a proper consideration of all subleading effects 
(vertex corrections, wave function renormalizations, the Meissner effect, etc.) would result in the cancellation of the 
gauge dependence in the order parameter, though it is hard to show this explicitly. Our consideration here just 
indicates that subleading corrections in the prefactor might be, though of order one, not small. 

The exponent and the power of gs in the prefactor in expression ( [42|) coincide with those in the expression for the 
gap obtained in Ref. ||20| by using the renormalization group method. Also, up to an overall factor of order one, this 
expression agrees with those in recent papers Refs. |p3,R3] where the SD equation in dense QCD has been studied by 
using somewhat different approaches. 

Regarding the validity of the expression in Eq. (p3), one has to bare in mind that most of the approximations 
made at the intermediate steps are adequate only for small values of the coupling constant. Therefore, in the most 
interesting case of the values of Ug ^ 0(1), the estimate in Eq. (E2) is not very reliable. Still we believe that it should 
give a correct order of magnitude for jA'-^'j. 

VI. CONCLUSION 

In this paper, we studied the phenomenon of color superconductivity in dense QCD with two light flavors in the 
framework of the SD equation in the improved rainbow approximation. Since our approach is based on the SD 
equations with a perturbative kernel, it may provide a reliable description of the system only at sufficiently large 
values of the chemical potential (say, larger than about a few GeV). Remarkably, the latter restriction on the range 
of validity is also sufficient for suppressing the non-perturbative contributions given by the instanton effects Q . 

Regarding our analytical dependence of the gap on the coupling constant presented in Eq. (B3) , a few comments 
are in order. First of all, because of the large prefactor in Eq. (|4^), the magnitude of the gap for realistic values of the 
chemical potential (200 to 700 MeV) and the coupling constant (0.8 to 0.4, respectively) can be as large as of order 
100 MeV (see Fig. p. Therefore, it is of the same order as the non-perturbative contribution due to the instantons. 



When the quark density is sufficiently large, however, the instanton effects become negligible and the order parameter 
is primarily determined by the gluon interaction. Indeed, the key point to notice is that the expression in Eq. (|42| ) 

contains as in the exponent in contrast to a^^ as appears in models with local four-fermion interactions induced, 
for example, by screened gluons or by instantons. Therefore, at small values of as, the estimate in Eq. (M3) is much 
larger compared to the estimates from four-fermion models. 

Thus, our analysis clearly shows that the long-range interaction, mediated by the unscreened magnetic modes of 
gluons, is responsible for the enhancement of the superconducting gap in the asymptotic region of high densities. This 
confirms the general conclusions of Ref. [ pO[ based primarily on renormalization group methods. 

It is natural to conjecture that expression (^), derived in the improved rainbow approximation, yields the exact 
essential singularity for the gap at as(/i)=0 (corresponding to the chemical potential /i going to infinity). In other 
words, the improved rainbow approximation can be a good leading approximation for large values of the chemical 
potential when the instanton contributions are suppressed. Indeed, the relevant region of momenta in the gap equation 
( p8| ) is |A(~)| ^ |q4| <^ fj,. Since as{fi) ~ l/ln(^/AQCD), expression (^) imphes that the gap |A(~)| goes to infinity 
with /i. Therefore the running coupling is small in that relevant region. This suggests that the exponent in expression 
(E3) for the gap is exact (on the other hand, the prcfactor can be influenced by higher order corrections in the kernel). 
This issue deserves further study. 

We note that the two approaches, i.e., using either the perturbative interaction or the non-perturbative instanton 
mediated one, seem to be complimentary ways of describing the intermediate region of dense QCD. While the stan- 
dard perturbative approach works well at large chemical potential, the instanton approach gives a more appropriate 
description at small chemical potential. 
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APPENDIX A: FORMULAE USED IN DERIVATION OF THE GAP EQUATION 



In derivation of Eqs. (p2|) - (|34|), we find it convenient to use the following representations for the kinetic terms of 
the quark propagators: 



{po + /i)7° + Bf=^° [{po -B\p\+ A.)A(+) + (po + B\p\+ A.)A(-) 
{po - m)7" + Bf^ 7" f(Po - B\p\ - A^)A(+) + (po + B\p\ - A^)A(-) 

[{po + t^h" + B^-' -7" 

Lf'O -r- J 

1 .,^^ 1 



I A(+) + i A(- 

Po + B\p\+f, P ^p,-B\p\+fi P 



[iPo-^ih" + Bfi]-' ^7" 



-Ai+) 



Ai-) 



Po + B\p\ - Ai ^ Po- B\p\ - A* ^ 



(Al) 
(A2) 

(A3) 
(A4) 



where the projectors Ap are defined in Eq. ( p5[ ) in the main text. 

In order to perform the angular integration in the right hand sides of the gap equation (p2|), one first need to 
calculate the following two types of traces over the Dirac indices: 



tr 
tr 



A(±)7'^A(±)7'' 



g'"'(l + i)-2<?^V°i + 



^.p^ + q^pf 



A(±)7'^Af )7''l = g'-'il -t) + 2g^^g''h - 



qf'p^ + 0^?^" 



(A5) 
(A6) 



where t — cos 9 is the cosine of the angle between three- vectors q and p, and irrelevant antisymmetric terms are 
denoted by ellipsis. 
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By contracting these traces with the projectors of the magnetic, electric and longitudinal types of gluon modes, we 
arrive at 



OW tr [a(±)7"A(^)7''' 



0(l) tr Al^h'A[^V 



O^ tr [a(±)7'^A(^)7''' 



0\f}, tr 



A(±)7'^A(±)7'' 



-2(l-t) 
= 2(l-t) 
= 2(l + t) 
= (l-t) 



q2 _^_ p2 _ 2qpt 

q^ + p^ + qpjl - t) 
g2 _j_ p2 _ 2qpt 

q^ + p'^ - qp{l + t) 
q2 ^ p2 _ 2qpt 

{q+pf + {qi-piY 



-(1-i) 



(1 + 



(g + p)' + (g4"P4)' 

q2 ^p2 _ 2qpt + (qr^ _ ^4)2 • 

(g - P)' + (g4 - P4)' 
q2 ^p2 _ 2qpt + (^4 - p4)2 ■ 



0|,t) tr [a(±)7^AW7''1 =(l + t) 



q2 _^_p2 _ 2qpt + (^4 - p4)2 ■ 

(g-p)^ + (g4 -Pif 

' q2 _^_p2 _ 2qpt + (^4 - p4)2 ■ 



(A7a) 
(A7b) 
(A7c) 
(A7d) 
(A7e) 
(A7f) 



where q = \q\, p = \p\, q^ = -iqo and pi = -ipo- 

Now, by making use of the explicit expressions ( A7), we easily perform the angular integrations of all types appearing 
in the gap equations: 



dn\q^p\ (1) N (±) u^(±) . 



^ , iq^-p' 

qp 



V3a'2-2„2 



■ arctan 



V3w/min(g,p) 



Lufq^p^ 



ujf + l?^ - P^l - w;max(g,p) 



(^q2 _p2-j2 -cjf ^^ a;; + |g + p| (g^ -p2)2-^f ^^ cjf + |g + pp - cj;|g + p| 



4 j^wf + |g + p|3 



Siufq'^p'^ uji + \q — p\ Quofq^p^ uof + \q — p\'^ — LOi\q — p\ 3qp uf + \q — p 



(A8) 



where ujf = {7r/2)M'^u) and uj = \qi — pi\. Similarly 
dn\q — p\ 



[- 

J w- 



■p\^+ujf 



OW tr A(±)7^A(T)7'' 



2 {q^^p^f+iof 

73^,2„2„2 



arctan 



V3ujiinm{q,p) 



qp 



icj fq^p^ 



Luf + \q'^ — p^l — w/max(g,p) 



(^2 _ p2)2 „ ^4 ^^ ^^ + \q+p\ ^ (g2_p2)2_^4 ^^ ^2 ^ |^ ^ ^(2 _ ^^j^ + p| 



4 ^^uf + \q+p\^ 



Sujfq'^p'^ u;i + \q — p\ dtofq^p'^ uf + \q — p\'^ — U}i\q — p\ 3qp ojf + \q — p\^ 



(A9) 



dn 



\q-p\^+uj^ + 2M^ 



-0(2) tr U(^^7^A(±)7''' 



27r 



qp 2g2p2 



(g2-p2)2 ^^ (g + p)2 
2M2 + Cj2 ^ (g - p)2 



[(g - pf + 2M2 + C^2] [(2M2 ^ ^2^2 ^ (^ ^ p)2^2] (g+p)2+2Af2 

2Af2(2Af2+w2) ^ 

[((? + pf + ^2] [(q - pf + cj2] (g + pf + u;2 
in 



{q - pf + 2Af 2 + t^2 



2M2 



(g — p)2 +0^2 



(AlO) 



/ 



dr^ 



9 - p12 + ^2 _^ 2M2 



0(1) tr [a(±'7'^A(^)7^' 



27r 



gp 2(72p2 



(^2 _p2)2 ^^ (g+p)2 



2M2 



(9 - P)' 



[{q + p)2 + 2M2 + t^2] [(2M2 + c^2)2 + (g_p)2^2] ^^ (^ ^ ^^2 ^ ^M^ ^ ^2 



2M2(2M2 + cj2) 
[{q + pf+u\{q-pf+uj^]^ {q+pf + 



{q - pf + 2Af 2 + tj2 



2A-/2 



In 



{q - pf + ^2 



(All) 



y ig-pi2- 



0(3) tr 



^2+c^2 P- 



A(±)7^A(±)7^ 



— nd 



2 
gp 



(g+p) 



2q2p2 



^ ' u^ [q+pf 



{q - pf + w2 



(A12) 
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dn 



|g-pT 



O'i^ tr [a(±)7^A(t)7^ 



— nd 



2 

qp 



{q — p)2 + cj2 



2g2p: 



2ti2 



(A13) 



As a result, the angular average of the gluon propagator (multiplied by q^ weight) in the vicinity of the Fermi 
surface is given by the following approximate expression: 



dD,V^^{q ~ p) tr 



A(±)7'^Af)7^ 



2i7r 



^m. 



(2^)= 



In- 



(2m)^ 



(e,-)2 + 2Af2+c^2 



(AM) 



where e~ = q — fx and d if the gauge fixing parameter. The three terms in the last expression are the leading 
contributions of the magnetic, electric and the longitudinal gluon modes, respectively. Because of the absence of the 
logarithmic factor in front of the gauge fixing parameter, the longitudinal gluon modes become relevant only in the 
next-to-leading order. 

Now, after performing the approximate integration over the spatial momentum, we arrive at 



q^dq 



4«7r2 



df7I?^,(g-p)tr[AW^^A(T)7" 



3\^r 



In ^-^ + - In ■ 



{2^,f 



2Af 2 + cj2 



3 ,\ 4i7r2 , A 

-d ~ zln— , 

2 J 3^ql + |A|2 c^' 



(A15) 



where A = e3''/2(2y^)V(V27rM5) = e^'^/H6{2TT)^/^iJ./{Nfas)^^^. 



APPENDIX B: THE SOLUTION OF THE GAP EQUATION 



In this Appendix we present a somewhat more rigorous solution of the integral equation ( p9| ) . 
In order to rewrite it in the form of a differential equation, we approximate the logarithm in the kernel by its 
asymptotes (below A = A^^' , p = p4^): 



A{p) 



2a. 



P dqA{q) , A , 2a, 
In ■ 



9^ 7o ^q^+W P 9^ Jp V?+^ 



dqA{q) A 

In—, 



(Bl) 



where, without loss of generality, A — e^'^/'^iA-n:)'^/'^ ^/ a^^"^ ^ appearing in the logarithm, is also used as an ultraviolet 
cut-off in integration over q. This equation, as is easy to check, is equivalent to the differential equation: 



pA"{j>) + A'(p) 



2a, 



A{p) 



97r ^p2 + |A|2 
along with the following infrared and ultraviolet boundary conditions: 

pA'{p)\^^,^Q (IR), A(A)^0 (UV). 



(B2) 



(B3) 



We solve the differential equation (B2) analytically in two regions p <^ \A\ and p ^ | A| and then match the solutions 
atp= |A|. 

In the infrared region p <C |A|. the solution that satisfies the IR boundary condition is given by the following Bessel 
function: 



(B4) 



A(p) = |A| Jo V 



where, by definition, v = ■y/SoT/QTr and the overall constant was chosen in such a way that A(0) = |A| 
In the other region, p^ |A|, the solution, consistent with the UV boundary condition, reads 

A{p) =Bsin ( -In— j . 

While matching the solutions and their derivatives at the point p = |A|, we get the relation: 



(B5) 



12 



Aexp 



and determine the value of the constant B: 



B^IA 



2 , /JoM 
- — arctan ^ , . 



v/w+^- 



(B6) 



(B7) 



The dependence (B6) of the order parameter on the chemical potential is presented in Fig. nl Here, to establish 
the function as(/x), we fix the magnitude of the QCD coupling constant by choosing as(1.76GeV) « 0.26. Then, by 
making use of the one- loop running of as , we obtain the functional dependence of as on /z. 

Finally, assuming that the coupling constant is small, i' ^ 1, we arrive at the analytical expression presented in 
Eq. (H): 



|A|~Aexp(-- + l + 0(i.2)) 



(B8) 
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FIG. 1. The dependence of the order parameter on the chemical potential. The function as(/i) is established by fixing 
as(1.76GeV) ^ 0.26 and using the one-loop running of the coupling constant. 
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